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I.  INTRODUCTION 


In  this  semi-annual  report  we  describe  progress  in  each  of  the  four  major 
areas  of  our  research.  These  areas  are  dynamic  routing,  multihop  packet  radio 
networks,  new  switching  techniques  in  networks,  and  algorithms. 

We  have  already  reported  on  our  proposed  dynamic  routing  technique.  In 
section  IIA  we  describe  the  results  of  simulations  that  have  confirmed  predictions 
about  performance  and  deepened  our  understanding  of  the  problem. 

Our  previous  work  on  multihop  packet  radio  has  been  extended  to  include 
the  effects  when  passive  acknowledgements  are  not  heard.  We  are  able  to 
extend  our  analytic  mociel  to  include  this  phenomenon  and  evaluate  the  perfor¬ 
mance  of  the  network  for  a  variety  of  protocols.  A  new  protocol  is  proposed 
that  essentially  guarantees  that  passive  acknowledgements  will  be  heard.  Details 
are  given  in  section  IIB. 

The  third  area  of  our  research  is  to  explore  new  switching  techniques  in 
networks.  There  are  many  applications  for  which  packet  switching  is  inefficient. 
Fortunately  modern  networks  have  the  resources  -  processing  power,  memory, 
distribution  of  intelligence  -  to  utilize  much  more  sophisticated  communication 
techniques.  We  have  been  exploring  this  subject  by  investigating  various 
techniques  and  applications  and  by  finding  bounds  on  the  best  that  can  be 
achieved.  In  section  IIC  we  describe  our  results  for  a  particular  application. 
Packet  switching  is  inefficient,  but  is  still  used  for  lack  of  an  alternative,  when 
a  large  number  of  users  transmit  short  messages  with  a  tight  time  delay  require¬ 
ment.  Here  we  find  an  optimum  packeting  strategy  and  compare  it  with  others. 
In  later  work  we  will  propose  alternative  switching  schemes  for  this  application. 

In  section  IID  we  describe  a  generalization  to  our  approach  to  algorithms. 
We  have  already  reported  on  many  algorithms  used  to  design  networks.  Here 
we  explore  the  asymptotic  and  probabilistic  behavior  of  algorithms. 

Research  in  each  of  these  four  areas  is  continuing. 
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II.  RESEARCH  SUMMARIES 

IIA.  A  Simulation  Study  of  A  Dynamic  Routing  Scheme 

1.  Introduction 

A  dynamic  routing  scheme  has  been  previously  proposed.  Instead 
of  assigning  fixed  paths  or  randomly  selecting  paths,  permissible 
paths  (e.g.,  minimum  hop  paths)  are  pre-assigned  to  each  node  pair. 
A  packet  arriving  at  a  node,  and  destined  for  another  node,  is 
assigned  to  one  of  two  types  of  queues.  A  dedicated  queue  is  a 
queue  whose  packets  are  served  by  one  and  only  one  dedicated  server. 
A  shared  queue  is  a  queue  whose  packets  may  be  served  by  two  or 
more  specified  servers.  A  server  is  equivalent  to  an  output  link. 
The  purpose  of  pre-selecting  paths  is  to  avoid  long  paths  and 
resultant  increased  network  traffic  that  would  occur  in  a  random 
routing  scheme.  The  purpose  of  the  shared  queue  is  to  let  a  portion 
of  the  nodal  traffic  experience  the  delay  of  a  multiple  server  queue 
and  therefore  reduce  the  average  packet  delay  at  a  node.  The  objec¬ 
tive  of  the  model  is  therefore  to  reduce  both  nodal  as  well  as  network 
average  packet  delay. 

We  have  conducted  a  series  of  simulation  studies  to  evaluate  the 
nodal  as  well  as  the  network  delay  under  the  proposed  dynamic 
routing  scheme.  Section  2  describes  the  nodal  and  network  models  as 
well  as  the  different  service  disciplines  used.  Section  3  gives  the  re¬ 
sults  of  nodal  and  network  delay  performance.  A  summary  is  given 
in  Section  4. 


2 .  The  Model 

The  Nodal  Model 

The  traffic  from  each  node  pair  is  preassigned  to  one  or  more 
pre-selected  paths  (and  therefore  server(s)  at  a  node).  Two  types 
of  queue  can  be  constructed  at  a  node.  A  dedicated  queue  can  only 
be  served  by  one  dedicated  server.  There  are  at  most  n  such  queues 
at  a  node  with  n  links.  A  shared  queue  can  be  served  by  any  of 
several  preassigned  servers.  There  could  be  as  many  as  2n-n-l 
shared  queues  at  a  node.  A  two  server,  three  queue  nodal  model  is 
shown  in  Figure  IIA.l. 

In  Figure  IIA.l,  traffic  arrives  at  the  different  queues  with 
arrival  rates  A .  (i  =  1,2,3)  packets/sec.  Queue  1  and  queue  3  are 
competing  with  server  1  for  service  while  queue  2  and  queue  3  are 
competing  for  service  by  server  2.  The  service  rate  is  p  packets/sec. 
If  one  of  the  competing  queues  (e.g.,  queue  1  or  queue  3)  is  empty 
(e.g.,  queue  1),  and  the  server  for  these  queues  is  idle,  a  packet  in 
the  non-empty  queue  is  served.  If  both  queues  are  not  empty  when 
the  server  has  just  completed  service,  a  pre-defined  service  discipline 
will  determine  which  queue  to  serve  next.  When  both  competing 
queues  are  empty  (e.g.,  queue  1  and  queue  3)  and  the  server  for 
these  queues  is  idle,  the  server  will  stay  idle  until  there  is  a  new 
arrival  to  one  of  these  queues.  The  idle  server  will  not  serve  the 
queue  which  is  not  pre-assigned  to  it  for  service.  This  is  the  major 
difference  between  our  nodal  model  and  a  multiple  server  queueing 
system.  The  operation  of  a  node  with  more  communications  links  is 
basically  the  same  as  that  of  a  two  link  node.  Shared  queues  will  be 
served  by  a  subset  of  the  servers. 
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FIGURE  IIA.l. 


A  NODAL  QUEUEING  MODEL  WITH  TWO  SERVERS 


Service  Disciplines 


The  nodal  model  described  in  the  previous  section  allows  different 
types  of  queues  to  compete  for  service.  A  service  discipline  must 
therefore  be  specified.  We  will  describe  four  different  service  disci¬ 
plines  for  the  nodal  model  of  Figure  IIA.l. 

a.  Priority  Service  Discipline  (PSD) 

A  server  always  searches  the  dedicated  queue  for  service 
prior  to  its  service  of  the  shared  queue.  For  example, 
server  1  in  Figure  IIA.l,  after  completing  a  service,  will 
always  look  to  queue  1  for  more  service.  Packets  in  queue 
3  will  be  served  by  server  1  only  if  queue  1  is  empty  when 
server  1  has  just  completed  a  service.  Therefore,  server  1 
is  busy  as  long  as  there  are  packets  to  be  served  in  queue 
1  or  queue  3.  Conversely,  server  1  is  idle  and  stays  idle 
when  both  queue  1  and  queue  3  are  empty.  The  instanta¬ 
neous  queue  length  of  queue  2  does  not  effect  any  of  the 
decision  rules  at  server  1  described  above.  A  packet 
entering  the  shared  queue  when  both  servers  are  idle  is 
served  randomly  by  either  server  1  or  server  2. 

b.  Alternate  Service  Discipline  (ASD) 

A  server  at  the  node  in  Figure  IIA.l,  whenever  possible, 
alternatively  serves  n  packets  from  the  dedicated  queue  and 
m  consecutive  ones  from  the  shared  queue.  The  server  is 
forced  to  serve  more  than  n  (or  m)  packets  from  the  dedi¬ 
cated  (or  the  shared)  queue  if  the  other  queue  happens  to 
be  empty.  Under  this  situation,  alternation  of  service  will 
begin  whenever  there  is  a  new  arrival  to  the  empty  queue. 
The  process  is  started  over  again  at  this  point. 
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c.  Serve  The  Longer  Queue  Discipline  (SLQD) 

The  queue  with  the  longer  queue  length  among  the  two 
competing  queues  (the  shared  and  the  dedicated  queues) 
gets  service  first.  A  server,  upon  a  service  completion, 
compares  the  queue  lengths  of  the  two  competing  queues, 
and  serves  a  packet  from  the  longer  queue.  The  SLQD 
tries  to  equalize  all  the  queue  lengths  at  a  node. 

d.  Random  Service  Discipline  (RSD) 

A  server  (e.g.,  server  1)  upon  completing  a  service, 
serves  packets  from  one  of  two  competing  queues  according 
to  the  following  algorithm:  it  serves  queue  1  if  queue  3  is 
empty  and  vise  versa;  it  randomly  selects  a  queue  if  both 
queue  1  and  queue  3  are  non-empty.  This  is  a  totally 
"uncontrolled"  service  discipline. 

The  Network  Model 

In  order  to  describe  the  dynamics  of  a  packet  switching  network 
employing  the  proposed  dynamic  routing  scheme,  one  must  (1)  con¬ 
struct  nodal  queueing  models  similar  to  the  one  in  Figure  IIA.l  for 
every  node  in  the  network,  (2)  pre-select  a  set  of  "permissible" 
paths  for  each  source-destination  node  pair,  and  (3)  assign  the  data 
flow  according  to  the  preselected  paths  between  adjacent  nodes.  In 
order  to  simplify  description  of  the  network  model,  we  have  chosen  a 
four  node  network  to  illustrate  the  dynamics  of  our  network  model. 

A  four  node  network  is  shown  in  Figure  IIA.2.  Every  node  is 
connected  to  two  links,  i.e.,  is  of  degree  two.  Therefore,  the  nodal 
queueing  model  at  each  node  is  exactly  the  same  as  that  of  the  two 
server,  three  queue  modei  described  previously.  The  permissible 


paths  between  sources  and  destinations  are  selected  based  on  a 
shortest  path  (minimum  hop)  algorithm.  Adjacent  nodes  will  only  use 
the  one  hop  path  between  them.  Non-adjacent  nodes  will  use  either 
one  of  the  two  two-hop  paths.  Therefore,  traffic  between  adjacent 
nodes  is  pre-assigned  to  the  dedicated  queues.  Traffic  between 
non-adjacent  nodes,  on  the  other  hand,  is  pre-assigned  to  the  shared 
queue  at  the  source  node  and  the  appropriate  dedicated  queue  at  the 
intermediate  node. 

Each  server  at  the  node  serves  the  two  competing  queues  (the 
dedicated  queue  and  the  shared  queue)  according  to  a  pre-specified 
service  discipline.  We  have  simulated  the  network  for  the  four  ser¬ 
vice  disciplines  discussed  previously.  A  service  discipline  where  the 
condition  of  the  next  node  is  "known"  to  all  servers  at  the  current 
node  and  used  to  determine  routing  is  also  studied. 
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3.  Results 

Simulation  programs  in  PL/1  have  been  written  to  evaluate  the 
delay  of  the  nodal  model  of  Figure  II A.  1  and  the  network  model  of 
Figure  IIA.2  operating  under  the  different  service  disciplines.  All 
simulation  results  are  based  on  statistics  collected  over  a  sufficient 
time  interval  so  as  to  represent  the  steady  state.  The  statistics 
collected  for  the  PSD  scheme  are  based  on  at  least  one  million  (some¬ 
times  four  million)  simulation  events  -including  both  arrivals  and 
departures.  In  order  to  conserve  computing  time,  the  statistics 
collected  for  other  service  disciplines  are  based  on  at  least  two  hundred 
thousand  simulation  events.  In  our  simulation  program,  all  arrival 
processes  to  source  nodes  are  assumed  to  be  Poisson  and  all  service 
distributions  are  assumed  to  be  exponential  and  independent  of  each 
other.  No  assumption  is  made  about  arrivals  processes  to  inter¬ 
mediate  nodes  in  the  network  model. 

Nodal  Delay  Performance 

The  percentage  of  the  traffic  that  is  shared  at  a  node  is  an 
important  parameter  affecting  the  nodal  delay,  We  have  studied  nodal 
delay  based  on  three  sets  of  nodal  traffic  rixes.  They  are  low 
(approximately  10%  to  15%),  moderate  (25%),  and  high  (50%)  percentages 
for  the  shared  traffic  at  a  node  respectively.  Fairly  extensive  simulation 
runs  were  conducted  throughout  the  whole  range  of  nodal  utilization 
(0  to  1)  for  the  PSD  and  a  moderate  percentage  of  shared  traffic. 

We  assumed  that  the  input  processes  to  the  node  were  Poisson. 
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TABLE  IIA.l 
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Nodal  Delay  Of  The  PSD  With  Low,  Moderate,  And  High 

Percentages  Of  Traffic  Being  Shared  (p-^1^2-^3  ,  p 
=  4  packets/sec) 


Shared  Arrival  Rate  to  Queue  Utilization  Nodal  Delay  (Seconds) 
(packets/sec. ) _ 


Traffic 

^•3 

P 

PSD 

M/M/2 

M/M/1 

2.8 

2.8 

1 

0.825 

1.01 

0.78 

1.43 

Low 

3.0 

3.0 

1 

0.875 

1.33 

1.07 

2.00 

3.2 

3.2 

1 

0.925 

2.08 

1.73 

3.33 

2.475 

2.475 

1.65 

0.825 

0.91 

0.78 

1.43 

Moderate 

2.625 

2.625 

1.75 

0.875 

1.18 

1.07 

2.00 

2.775  2.775 


1.85  0.925  1.87  1.73  3.33 
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FIGURE  IIA.3.  DELAY-THROUPUT  RELATIONSHIP  OF  NODAL  MODEL 
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Simulation  Results  Of  The  Priority  Service  Discipline 

The  effects  of  the  percentage  of  the  traffic  that  is  shared  on 
nodal  delay  are  tabulated  in  Table  IIA.l.  Notice  that  the  higher  the 
percentage,  the  closer  the  nodal  delay  is  to  that  of  an  M/M/2  queueing 
system.  Figure  IIA.3  plots  the  delay-throughput  relationship  of  the 
PSD  for  different  percentages  of  shared  traffic.  The  delay  for  the 
M/M/2  and  the  M/M/1  queueing  systems  are  also  plotted  for  reference. 

In  Figure  IIA.4,  a  measure  of  the  "closeness"  of  the  PSD  to 
M/M/2  queueing  system  performance,  is  plotted  against  the  percentage 
of  shared  traffic  for  a  node  operated  under  the  PSD  scheme.  Two 
utilizations  are  studied.  The  ordinate  is  the  percentage  of  multiple 
server  performance  that  is  obtained  by  the  PSD.  At  high  traffic 
utilization  (p=0.875),  only  25%  of  the  nodal  traffic  need  to  be  shared 
to  achieve  90%.  At  higher  nodal  utilization,  the  same  effect  is 
reached  with  less  shared  traffic.  At  moderate  nodal  traffic  utilization 
(p=0.5),  th  PSD  achieves  65%  of  the  multiple  server  performance  with 
only  25%  of  the  traffic  being  shared. 

In  Figure  IIA.5,  we  have  redrawn  the  plot  of  Figure  IIA.3  for 
the  average  nodal  delay  in  the  range  from  0  to  0.5  seconds  for  the 
case  of  25%  shared  traffic.  For  a  fixed  delay  of  0.5  second,  the 
utilization  of  the  PSD  scheme  is  about  25%  greater  than  that  of  the 
M/M/1  model  (a  commonly  used  fixed  routing  nodal  nodel).  Thus  for 
the  same  throughput  a  significant  decrease  in  required  channal 
capacity  is  obtained. 

Results  For  Other  Service  Disciplines 

In  this  section,  we  discuss  nodal  delay  for  other  service  disci¬ 
plines.  For  these  disciplines,  we  only  studied  the  nodal  performance 
for  p=0.875  and  50%  of  the  traffic  being  shared. 


Table  IIA.2  tabulates  the  nodal  delay  from  our  nodal  simulation  for 

the  ASD,  the  SLQD,  and  the  RSD.  We  can  make  the  following  obser¬ 
vations  from  Table  IIA.2: 

a.  A  general  characteristic  of  the  results  in  Table  IIA.2  is  that:  at 
high  utilizatons  (p=0.875  in  this  case)  with  a  high  percentage  of 
the  traffic  shared  (50%),  the  nodal  delay  is  closer  to  that  of  the 
M/M/2  than  the  M/M/1  queueing  system  regardless  of  the  service 
disciplines . 

b.  The  PSD  has  the  best  nodal  delay  among  all  service  disciplines. 

c.  The  nodal  delay  of  the  SLQD,  with  the  modification  of  not  serv¬ 
ing  (whenever  possible)  two  packets  in  a  row  from  the  shared 
queue,  is  the  closest  to  that  of  the  PSD. 

d.  The  ASD  with  parameter  n  =  1  and  the  threshold  at  the  shared 
queue  set  to  zero,  has  almost  the  worst  nodal  delay  performance. 

e.  The  RSD,  being  that  it  is  an  "uncontrolled"  service  discipline, 
does  perform  the  worst  among  all  service  disciplines.  However, 
the  RSD  performs  better  than  the  M/M/1  queuing  system  by 
almost  40%. 
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TABLE  IIA.2  Nodal  Simulation  Results 

(\1=A2=l-75  packets/sec,  A3=3.5  packets/sec, 
(j=4  packets/sec,  and  p=0.875) 


Service  Discipline 


PSD 


Nodal  Delay 
(Seconds) 


Serve  the  Longer 
Queue  but  not  y 
consecutive  "shared 
ueue"  packets 


SD 

Alternate  Service: 
one  from  shared 
queue  and  n 
from  dedicated 
queue  but  do  not 
serve  shared 
queue  when  its 
queue  length  <  TH 


RSD  Random 


M/M/2 


_  M/M/1 


1.09 


1.08 


1.13 


1.11 


=  2 

1.12 

=  3 

1.05 

n  =  2 


=  3 


n  =  2 


=  3 


1.09 


1.05 


1.07 


1.08 


1.05 


1.04 


1.06 


1.07 


1.05 


1.13 


1.01 


2.00 
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Network  Delay  Performance 

Analysis 

In  this  section,  we  study  the  average  network  delay  performance. 
In  our  routing  strategy,  instead  of  treating  each  communications  link 
separately,  we  focus  our  interest  on  the  whole  node,  including  all  its 
communications  links.  As  a  result,  our  focus  is  on  nodal  delay  as 
opposed  to  link  delay,  which  is  usually  studied. 

Let  us  assume  that  we  have  an  expression  for  nodal  delay  at 
node  j,  which  is  denoted  by  Tj.  The  average  number  of  packets  N. 
at  node  j  can  therefore  be  written  as  Nj  =  NjTj,  where  A!,  is  the  total 
traffic  rate  (from  both  external  and  internal  traffic  sources)  at  node 
j.  The  average  number  of  packets  in  the  network  (N)  is  therefore 
the  sum  of  the  average  number  of  packets  at  each  node.  That  is, 
N  =  all  fiodes  Nj'  The  avera9e  network  delay  (T)  is  therefore 

T  *  \  alfnodes  <IIIW> 

where  Y  =  j  ^  Yj^  is  the  total  external  network  traffic  demand  and 
is  the  offered  rate  of  traffic  between  pairs  (j,k).  Equation  (IIA.l)  is 
therefore  a  general  expression  for  the  average  network  delay. 

To  analyze  the  network  delay  we  make  the  following  assumptions. 
We  assume  the  node  models  used  previously.  The  traffic  is  found  by 
applying  the  routing  implied  by  the  preassigned  paths.  We  can  then 
calculate  the  rate  in  packets  per  second  of  arrivals  to  the  shared  and 
dedicated  queues.  In  our  network  analysis  we  assume  these  arrivals 
to  be  a  Poisson  process  and  compute  the  nodal  delay  for  each  node 
using  the  techniques  described  above.  Equation  (IIA.l)  gives  us  the 
network  delay.  Below  we  compare  this  technique  with  the  results  of 
two  simulations.  The  first  dispenses  with  the  Poisson  assumption  and 
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accurately  models  arrivals  to  intermediate  nodes  in  paths.  To  make 
the  simulation  simpler,  packet  lengths  are  independently  reassigned  at 
each  node  in  a  path.  The  second  simulation  dispenses  with  this 
approximation  and  retains  the  initial  packet  length  for  a  packet 
throughout  its  path. 

Simulation 

In  this  section,  we  will  focus  our  discussion  on  two  PL/1  simula¬ 
tions  which  were  written  to  evaluate  the  network  delay  performance. 
The  interarrival  time  of  external  arrivals  to  network  nodes  is  assumed 
to  be  exponentially  distributed.  The  packet  length  is  assumed  to  be 
exponentially  distributed.  All  servers  are  assumed  to  have  the  same 
service  capacity  (for  convenience  chosen  as  unity).  Therefore,  the 
average  service  rate  of  each  server  is  the  same.  In  most  of  our 
simulation  studies,  we  have  used  p=0.875  ,  50%  shared  traffic  at  the 
nodes,  and  a  symmetrical  network.  Other  traffic  utilizations,  nodal 
traffic  mixes,  and  asymmetric  network  traffic  demands  are  also  studied 
but  not  as  extensively  as  the  above. 

The  first  simulation  allows  the  packet  lengths  to  be  reassigned  at 
the  next  node.  It  does  not  make  any  assumptions  about  the  output 
process  from  a  server.  A  server  serves  its  competing  queues  accord¬ 
ing  to  a  pre-specified  service  discipline.  After  the  server  finishes 
service  of  a  packet,  the  packet  either  exits  from  the  network  or  joins 
the  next  queue  for  more  service.  Once  it  arrives  at  the  next  queue 
(in  the  next  node),  a  new  packet  length  from  the  same  exponential 
distribution  is  assigned  to  it.  It  then  is  put  on  the  appropriate 
queue  and  waits  for  service.  We  compare  the  results  of  our  network 
simulation  with  our  analysis  based  upon  nodal  simulation  results.  In 
the  latter  we  assume  all  inputs  to  the  node  to  be  Poisson. 
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The  second  simulator  simulates  the  true  service  situation  in  an 
operating  packet  network.  That  is,  a  particular  packet  when  served 
by  the  next  server  will  require  the  same  amount  of  service  time. 
When  the  results  obtained  from  the  second  simulation  are  compared 
with  our  analysis  by  using  our  nodal  simulation  results,  we  can 
evaluate  the  effect  of  both  the  Poisson  output  stream  and  the  packet 
length  reassignment  assumptions. 

The  First  Network  Simulation 
Local  Service  Disciplines 

All  four  service  disciplines  are  studied.  The  simulation  results 
in  Table  IIA.3  and  in  previous  work  reveal  the  following: 

a.  The  average  network  delay  from  the  network  simulation  is  signif¬ 
icantly  colser  to  that  of  the  M/M/2  nodal  model  than  to  the 
M/M/1  nodal  model. 

b .  The  difference  between  the  analytical  results  using  the  nodal 
model  and  those  obtained  from  the  network  simulation,  for  the 
PSD,  increases  as  the  percentage  of  shared  traffic  increases  at  a 
node. 

c.  With  the  same  traffic  mix  at  a  node,  the  difference  between  the 
network  delay  found  analytically  and  those  obtained  from  the 
network  simulation  for  the  PSD  increases  as  the  traffic  load 
increases. 

d.  The  network  delay  obtained  from  the  network  simulation  for  the 
Alternate  Service  Discipline  (ASD)  and  the  Serve  the  Longer 
Queue  Discipline  (SLQD)  are  lower  than  that  of  the  PSD  scheme. 
The  difference  in  these  delay  are  at  most  11%.  The  network 
delay  of  the  Random  Service  Discipline  (RSD)  scheme  is  approxi¬ 
mately  the  same  as  that  of  the  PSD  scheme. 


e.  The  difference  between  the  results  obtained  analytically  and  from 
our  network  simulation  for  other  service  disciplines  is  substan¬ 
tially  large  (12  to  26  percent  difference). 

f.  The  parameter  n  (the  number  of  packets  to  be  served  from  the 
dedicated  queue  before  the  server  alternates  its  service)  in  the 
ASD  has  an  effect  on  network  delay. 

g.  In  the  ASD,  the  threshold  parameter  set  at  the  shared  queue  to 
avoid  the  shared  queue  becoming  empty  before  the  dedicated 
queues,  also  has  an  effect  on  network  delay. 

h.  The  network  delay  of  the  RSD  -  an  "uncontrolled"  and  "unpre¬ 
dictable"  service  discipline  -  is  about  the  same  as  that  of  the 
PSD  and  the  "worst  cases"  of  the  ASD  and  the  SLQD. 

Although  we  are  operating  at  high  utilization,  we  are  getting 

most,  but  not  all,  of  the  expected  improvement.  Many  different 
nodal  strategies  have  only  a  modest  effect  on  resolving  this  gap.  We 
find  that  it  was  easy  to  generate  strategies  that  produced  multiple 
server  behavior.  But  some  of  them  have  a  profound  effect  on  the 
output  process  from  the  server,  and  thus  affect  the  Poisson  assump¬ 
tion  for  inputs  to  the  next  node.  The  PSD  is  notorious  in  this  regard 
Usually  any  traffic  one  hop  away  from  its  destination  will  be  placed  in 
a  dedicated  queue.  Thus,  only  continuing  traffic  will  be  found  in 
nondedicated  queues.  These  have  lower  priority  and  will  be  trans¬ 
mitted  only  when  the  dedicated  queue  is  empty.  But  then  a  burst  of 
several  packets  is  likely  to  be  sent.  A  bursty  arrival  process  results 
in  much  poorer  queue  performance  than  a  Poisson  process.  This 
contributes  significantly  to  the  gap  in  delay.  The  design  of  a  nodal 
strategy  must  consider  not  only  multiserver  performance  at  a  node, 
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• 

TABLE  IIA.3 

Network  Simulation 
(Xil=Ai2=0  packet/sec,  A.^ 
(jil=|jp=4.0  packets/sec. 

Results 

3.5  packets/sec, 

3=0.875) 

• 

Network  Delay  (seconds) 

Service  Discipline 

Nodal  Model 

Simulation 

%  Difference 

PSD 

2.06 

2.75 

33% 

y 

n  =  1 

2.22 

2.65 

20% 

ASD 

TH  =  0  n  =  2 

2.24 

2.51 

12% 

- 

Alternate  Service: 

n  =  3 

2.10 

2.58 

23% 

•; 

one  from  shared 

n  =  1 

2.18 

2.50 

15% 

queue  and  n 

TH  =  1  n  =  2 

2.10 

2.50 

19% 

from  dedicated 

n  =  3 

2.14 

2.53 

18% 

queue  but  do 

n  =  1 

2.16 

2.59 

20% 

*•1 

not  serve 

TH  =  2  n  =  2 

2.10 

2.52 

20% 

shared  queue 

n  =  3 

2.08 

2.62 

26% 

when  its  queue 

n  =  1 

2.12 

2.51 

18% 

length  <  TH 

TH  =  3  n  =  2 

2.14 

2.49 

17% 

'•I 

n  =  3 

2.10 

2.55 

21% 

SLQD 

y  =  1 

2.08 

2.45 

18% 

Serve  the  Longer 

V  =  2 

2.18 

2.58 

19% 

Queue  but  not  y 

y  =  3 

2.16 

2.67 

24% 

..  -m 

consecutive 

y  =  co 

2.26 

2.67 

19% 

-  ■  -1 

shared  queue 
packets 

RSD 

2.26 

2.71 

20% 

-  ■~®4 

M/M/2 

2.02 

- 

- 

. 

M/M/1 

4.00 

- 

- 

: 

* 

%  Difference  = 

T  NETWORK  (Simulation)  ’  T 

NETWORK  (Nodal  Model) 

: 

T 

NETWORK  Model) 

.  "»4 
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but  the  resultant  output  process  as  well.  The  usual  way  such  a 
burst  could  occur  is  when  the  packets  are  short.  But  these  should 
have  made  the  next  node's  job  easier.  'Jnfortunately ,  the  second 
aspect  of  the  independence  assumption  now  comes  into  play  -  packet 
lengths  are  reassigned  at  each  node.  This  accentuates  the  effect  of 
burstiness.  When  the  simulation  is  altered  to  take  into  account  de¬ 
pendence  between  nodes,  hence  dispensing  with  the  independence 
assumption,  the  difficulty  disappeared. 

A  Look  Ahead  Service  Scheme 

A  common  situaton  of  all  the  "local"  service  disciplines  (The 
PSD,  the  ASD,  the  SLQD,  and  the  RSD)  is  that  a  packet  may  be 
waiting  to  be  served  in  the  shared  queue  at  one  node  while  there  is 
an  idle  server  at  the  next  node.  Conversely,  a  packet  may  wait  in 
the  dedicated  queue  for  service  when  there  is  no  need  to  serve  a 
shared  queue  packet  because  the  next  server  has  plenty  of  packets  to 
serve  at  the  next  node.  Both  phenomena  cause  degradation  in  net¬ 
work  performance.  Both  are  due  to  the  fact  that  each  server,  upon 
completing  a  service,  does  not  know  what  the  next  server  wants. 
Consequently,  the  server  may  serve  a  "wrong"  packet  from  time  to 
time.  Therefore,  it  is  not  the  selection  of  a  queue  but  the  timing  of 
this  selection  process  which  influences  the  service  decision.  In  all  of 
the  "local"  service  disciplines  which  we  have  discussed  so  far,  the 
timing  factor  was  not  considered.  We  have  designed  a  look  ahead 
service  discipline  and  have  examined  its  performance.  The  basic 
algorithm  is  shown  in  Figure  II A. 6. 


Table  IIA.4  shows  our  results.  We  make  the  following  observations: 

a.  The  network  delay  of  the  best  look  ahead  scheme  has  an  improve¬ 
ment  of  about  11%  over  that  of  the  best  local  scheme. 

b.  The  look  ahead  service  scheme  performs  consistently  better  than 
the  local  service  disciplines. 

c.  When  both  servers  are  empty,  it  is  better  to  serve  an  arrival  to 
the  shared  queue  on  the  basis  of  downstream  information,  than 
to  do  it  randomly. 

The  Second  Network  Simulation  -  No  Packet  Length  Reassignment 

The  analytically  obtained  network  and  our  first  simulation  differ 
not  only  in  the  PSD  but  also  in  the  ASD  and  the  SLQD  where  the 
output  stream  should  be  somewhat  smoother.  Furthermore,  the  packet 
iength  reassignment  assumption  does  not  work  well  in  a  small  network 
like  our  four  node  network  with  very  few  traffic  streams  flowing  at 
each  node.  This  motivated  us  to  build  the  second  network  simulator 
whereby  the  packet  length  reassignment  assumption  is  not  made. 


(a)  A  PART  OF  THE  FOUR  NODE  NETWORK 


X2  ■  Qa  OR  Q2i+  i/2  Q23 

(b)  TVE  LOOK  AHEAD  ALGORITHM 

FIGURE  IIA.6  A  LOOK  AHEAD  SERVICE  SCHEME 
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TABLE  IIA.4  Network  Simulation  Results  Of  A  Look  Ahead  Routing  Strategy 
(.Vl=Ap=0  packet/sec,  A. 2=3. 5  packets/sec,  p. ^=jj. 0=4  packets/sec, 
p=0.875) 


Threshold 


0  (PSD) 


L 


2 


3 


4 


5 


6 


Network  Delay  (seconds) 

X2  =  Q21*  X2  =  Q21**  X2  =  Q21  +  ^23** 
2.75  2.75 


2.33  2.26  2.27 


2.36  2.29  2.19 


2.37  2.32  2.26 


2.35  2.42  2.32 


2.35 


2.41 


When  both  Sn  and  S12  are  not  busy  and  there  is  an  arrival  to  Qx3,  this 
packet  will  be  served  by  either  S1]L  or  S12. 

Under  the  same  situation  as  above,  X2  will  tie  compared  with  X4  (the  same 
parameter  at  node  4). 

The  packet  in  Q  13  will  be  served  by  SX1  (or  S12  )  if  X2  <  X4  (or  X2  >  X4) 


Dynamic  Routing  Scheme  With  Locai  Node  Service  Disciplines 

We  have  used  our  second  simulation  to  evaluate  the  network 
delay  performance  under  the  various  service  disciplines  (PSD,  ASD, 
SLQD,  and  RSD)  discussed  previously.  Each  node,  operates  indepen¬ 
dently,  serves  its  queues  according  to  the  pre-specified  service 
discipline.  Traffic  information  at  each  node  is  not  exchanged  between 
adjacent  nodes.  Therefore,  only  "local"  node  service  disciplines  are 
discussed  in  this  section.  We  have  also  examined  the  average  network 
delay  under  two  fixed  routing  strategies. 

Table  IIA.5  tabulates  our  results.  We  make  the  following  observations: 

a.  Removing  the  packet  length  reassignment  assumption  greatly 
improves  the  accuracy  of  our  analytic  nodal  method. 

b.  The  packet  length  reassignment  assumption  has  more  effect  on 
the  network  delay  performance  for  the  PSD  than  any  other 
service  disciplines. 

c.  The  PSD  is  the  best  "local"  service  disciplines  among  all  other 
local  service  disciplines. 

d.  The  network  delay  of  our  dynamic  routing  scheme  under  any  one 
of  the  service  disciplines  is  significantly  better  than  either  of 
the  benchmark  fixed  routing  strategies. 

Dynamic  Routing  With  Look  Ahead  Service  Scheme 

The  "clock  driven"  network  simulation  keeps  track  of  how  much 
service  time  a  packet  requires  from  a  server.  This  enables  a  server 
to  know  how  much  work  (in  terms  of  service  time)  the  next  server 
has  at  any  given  instant  of  time.  A  server  can  make  an  accurate 
decision  to  keep  the  network  from  losing  its  service  capacity.  A 
server,  upon  completing  a  service,  compares  the  total  amount  of  work 
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TABLE  IIA.5  Network  Simulation  Results  For  Different  Service  Disciplines  - 
Local  Decision  Rules 

(Ail=Ai2=0  packet/second,  A*2=3 . 5  packets/second,  Mj^=Mp=4 
packets/second,  p=0,875) 


Network  Delay  (seconds) 

Service 

First 

Second 

Analytic 

Discipline 

Simulation 

Simulation 

PSD 

2.75 

2.16 

2.06 

ASD 

2.65 

2.38 

2.22 

2.45 

2.32 

2.08 

RSD 

2.73 

2.56 

2.26 

Fixed  Routing 
#1 

-- 

3.84 

— 

Fixed  Routing 
#2 

-- 

3.26 

-- 

at  the  next  server  (the  remaining 

service  time 

for  the  packet  still  in 

service  at  the  next  server,  the 

total  service 

time  in  the 

dedicated 

queue  at  the  next  node,  the  total  service  time  in  the  shared  queue  at 
the  next  node,  and  a  threshold  parameter  which  incorporates  the  fact 
that  the  shared  queue  packets  are  served  by  both  servers  at  the 
next  node  and  there  are  new  arrivals  to  the  next  node  during  the 
service  time)  to  the  sum  of  the  service  times  of  the  first  packet  in 
each  of  its  competing  (shared  and  dedicated)  queues.  The  server 


II  A.  27 


serves  the  shared  queue  packet  if  the  total  amount  of  work  at  the 
next  server  is  less  than  the  sum  of  the  service  times  of  the  first 
packet  in  each  of  its  competing  queue.  Should  the  server  not  serve 
the  shared  queue  packet  immediately,  the  next  server  could  be  idle 
for  a  period  of  time  equal  to  the  difference  between  these  two  quan¬ 
tities.  The  server  would  otherwise  serve  the  dedicated  queue  packet. 
This  look  ahead  service  scheme  would  be  ideal  in  our  dynamic  routing 
if  the  threshold  parameter  could  be  correctly  set.  The  network  delay 
with  both  exponentially  distributed  and  constant  packet  lengths  are 
evaluated.  The  following  observations  can  be  made  (see  Table  IIA.6): 

a.  The  PSD  is  equivalent  to  setting  the  threshold  value  to  negative 
infinity. 

b.  The  network  delay  is  insensitive  to  the  threshold  values  near  its 
"optimal"  setting . 

c.  The  difference  between  the  network  delay  of  the  "best"  local 
service  discipline  (the  PSD)  and  the  look  ahead  service  scheme 
with  the  "best"  choice  of  the  threshold  value  is  about  10%. 

Assymetric  Traffic  Demands  And  Service  Rates 

We  have  illustrated  significant  improvements  in  network  delay 
using  our  adaptive  techniques  for  the  four  node  network.  However, 
only  symmetric  traffic  loads  have  been  studied  so  far.  In  order  to 
demonstrate  the  robustness  of  our  routing  strategy,  a  certain  degree 
of  asymmetry  was  introduced  into  the  network.  We  have  evaluated 
the  network  delay  of  the  four  node  symmetric  network  with  asym¬ 
metric  traffic  and  with  different  link  capacities. 

We  make  the  following  observations  from  the  results  of  Tables 


IIA.7  and  II A. 8: 


The  RSD  performs  the  worst  among  all  service  disciplines;  the 
PSD  performs  the  best  among  all  routing  strategies  based  on 
local  traffic  information;  the  "best"  look  ahead  scheme  performs 


TABLE  IIA.6  Network  Simulation  Results  For  A  Look  Ahead  Routing  Strategy 
(\il=\i2=0  packet/second,  Aj2=3.5  packets/second, 

|Jil=Pi2=4  packets/second,  p=0/875) 


Network  Delay  (seconds) 


Local 


Routing  Strategy 


,  PSD 


W 


Exponential  Packet 
Lengths 


Constant  Packet 
Lengths 


2.56 

1.31 

.25 

1.9 

6 

1.17 

2.13 
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TABLE  IIA.7  Network  Delay  With  Asymmetric  Traffic  Demands 


(A. -j=2  packets/sec,  A^l  Packet/sec,  ^3=2  packets/sec, 
p-=4  packets/sec,  p=0/875) 


Routing  Strategy 

Network  Delay  (seconds) 

PSD 

1.91 

Local  SLQD 

1.96 

RSD 


2.67 
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TABLE  IIA.8  Network  Delay  With  Asymmetric  Service  Rates 


(Xil=Xi2=0  packets/sec,  X.^3.5  packets/sec, 
M 1 i=M22=M3i=M42=5  packets/sec, 
Mi2=M2i=iJ32=M4i=3  packets/sec, p=0. 875) 


Routing  Strategy 

Network  Delay  (seconds) 

Local 

PSD 

2.23 

-0.25 

2.07 

Look 

0 

2.03 

TH 

0.25 

2.04 

Ahead 

0.75 

2.15 

1.00 

2.22 

better  than  the  PSD  in  terms  of  average  network  delay;  and  the 
threshold  value  set  in  the  look  ahead  scheme  is  insensitive  to  the 
network  delay  peformance  over  a  wide  range. 

4.  Conclusion 

We  have  studied  a  dynamic  routing  scheme  for  packet  switching 
networks.  Its  network  delay  performance  is  superior  to  that  of  fixed 
routing.  We  have  also  developed  and  evaluated  an  analytic  technique 
for  evaluating  adaptive  routing  performance.  Extensive  simulations 
confirm  these  results.  A  locally  adaptive  scheme  performs  almost  as 
well  as  one  that  has  more  information. 
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The  Effect  of  Imperfect  Acknowledgments  on  the 
Throughput  of  Packet  Radio  Networks. 


1 .  Introduction 

In  this  section  we  extend  our  previous  work  on  throughput  in 
multihop  packet  radio  networks  to  include  the  effects  of  passive 
acknowledgments  not  being  heard.  We  also  evaluate  the  performance 
of  a  new  protocol  which  greatly  increases  the  probability  of  hearing 
these  acknowledgments. 

All  terminals  (nodes)  operate  under  a  CSMA  protocol.  Node  i 

will  transmit  a  scheduled  packet  to  j,  a  neighbor  of  i,  if  i  and  all  i's 

neighbors  are  idle.  Denote  the  set  of  all  neighbors  of  i,  including  i 

itself,  as  N..  The  transmission  will  be  heard  by  j,  assuming  perfect 

capture,  if  all  nodes  in  Nj  are  also  idle.  Let  P(A)  be  the  probability 

that  all  nodes  in  the  set  A  are  idle.  Let  G.  .  be  the  rate  of  sched- 

*  f  J 

uled  traffic  from  i  to  j.  Let  S.  •  be  the  rate  of  successfully  trans- 

*  /  J 

mitted  packets  from  i  to  j.  Then  since  the  probability  of  transmission 
is  P(Nj),  and  the  probability  of  successful  reception  is  P(Nj|N.),  we 
have 


Si,j  =  Gi/<rVNi>P<Ni>  =  Gi,jP(1VV 


(IIB-1) 


When  j  transmits  the  packet  to  the  next  node  in  the  route,  i  can 
hear  the  transmission,  and  thus  receive  a  passive  acknowledgment,  if 
all  of  its  neighbors  (except  j)  are  idle.  If  it  does  not  hear  the 
transmission,  then  after  a  suitable  time-out  period  it  retransmits  the 
packet.  Previously  we  assumed  all  passive  acknowledgments  were 
heard.  Here  we  see  that  a  passive  acknowledgment  is  heard  with 
probability  P(N.|Nj).  We  assume  that  the  time-out  period  is  selected 
so  that  all  nodes  have  the  opportunity  of  hearing  exactly  one  trans- 


mission.  If  they  fail  to  hear  it,  the  packet  is  retransmitted.  This  is 
repeated  until  an  acknowledgment  is  heard.  The  assumption  that 
exactly  one  transmission  may  be  heard  in  a  time-out  period  is  equiva¬ 
lent  to  assuming  that  the  same  time  out  period  is  used  for  all  nodes 
and  that  it  is  long  enough  to  insure  one  transmission. 

The  rate  S-  ■  in  equation  (IIB-1)  includes  repeated  packets, 

*  •  J 

because  of  failure  to  hear  acknowledgments.  Let  Sf  •  be  the  rate  of 
packets  for  which  passive  acknowledgments  are  heard.  Then 

sf,i  ■  <i.j  si,j  (I1B-2) 

P(N,,Nj) 

where  q.j  =  P(Nj|Nj)  =  '  (IIB-3) 

We  have  already  shown  how  to  solve  sets  of  equations  like 

(IIB-1  ^  for  the  G.  .  as  a  function  of  the  S-  .,  and  found  the  maximum 

*  /  J  w] 

possible  S.  .  which  we  defined  as  the  throughput.  For  example,  in  a 

chain,  all  S.  .  =  S  for  j  =  i  and  i  -  1.  Here  we  see  how  to  modify 
*  /  J 

this  procedure  to  account  for  passive  acknowledgments .  The  S?  j  are 
determined  from  the  offered  traffic  and  the  routing. 

We  identify  two  situations.  Packets  are  retransmitted  if  they  are 
unsuccessfully  transmitted  or  if  the  passive  acknowledgment  is  not 
heard.  In  the  latter  case,  a  duplicate  packet  is  sent.  If  the 
receiving  node  recognizes  the  duplicate  packet  and  deletes  it,  then 
Sjf  j  is  the  desired  rate  of  traffic  between  i  and  j.  If  a  simpler  node 
does  not  have  this  capability,  then  duplicate  packets  are  retrans¬ 
mitted,  and  after  several  hops  many  duplicates  exist.  In  this  latter 
case,  Sf  .  includes  these  duplicate  packets  and  must  be  related  to  the 
desired  rate  of  traffic.  We  will  first  consider  a  chain. 
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2.  A  chain  network 

Consider  the  chain  shown  in  figure  IIB-1  consisting  of  n  nodes. 
Assume  that  duplicate  packets  are  not  detected. 


0 - 0 - 0— . — 0 - 0 

123  n-1  n 


Figure  IIB-1.  A  chain  network 


We  assume  that  node  1  and  n  wish  to  send  S  packets  per  second  to 
each  other.  All  other  nodes  are  just  repeaters.  When  node  2  relays 
a  packet  to  node  3,  node  1  must  be  idle,  and  all  passive  acknowledg¬ 
ments  are  always  heard.  Thus  2  =  ^  anc*' 

S*,2  *  S1.2  =  S 

At  the  other  end,  node  n  does  not  repeat  the  packets,  so  passive 
acknowledgments  are  not  possible.  However  every  transmission  must 
ba  successful  since  there  are  no  interfering  nodes.  Thus 


S£-l,n  =  Vl,n 


(IIB-5) 


Node  i  attempts  to  transmit  all  S.  ,  .  successful  packets  to  i+1. 

I™ X  ,  1 

Thus 


V.iH  =  Si-u  '  1  =  2 . 

Similarly,  in  the  reverse  direction,  we  have 


S*  =  S  1  =  S 
n,n-l  n,n-l 

^2,1  =  ^2,1 

and  SM.i  =  Si,i-1  • 


i  =  2, . . . ,n-l 


(IIB-6) 


(IIB-7) 


r- 

"  •  '  • . . ' 

• 

* 
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« 

• 

Equations  (IIB-4  to  7)  relate  Sj  j  and  S;  ■*  to  the  desired  traffic,  S. 

fa 

Note,  by  symmetry,  we  obtain  relations  like  equations  (IIB-7)  by 

• 

replacing  i  and  j  by  n-i+1  and  n-j+1,  respectively. 

• 

As  an  example  we  evaluate  the  throughput  of  a  four  node  chain. 

m 

First  we  assume  perfect  reception  of  acknowledgments  and 

use  only 

• 

equation  (II3-1)  with  all  S-  .=S.  Then 

*  f  J 

S  =  G12P(1,2,3)  =  G23P(1. 2,3,4)  =  G21P(1,2,3) 

(IIB-8) 

We  let  Gl2  =  Gj,  G21  +  G23  =  G2,  and  note  by  symmetry  that 

G3  —  G2 

• 

• 

G4  =  G1(  and  P(l,2,3)  =  P(2,3,4).  From  our  previous  work 

we  have 

> 

1+G4  1+Gx 

P(l,2,3)  =  — p-'= 

-#! 

P(l,2,3,4)  =  i 

(IIB-9) 

- 

n 

D  =  1+2G1+2G2+2G1G2+Gt  =  (l+Gl)%2G2(l+G1) 

■•j 

Thus  we  get 

•  * 

m 

G2  =  G,.  +  GiCl+GO  =  G1(2+G1) 

-•< 

j; 

Gx  Gx 

1+G!+2G2  (1+G1)+2G1(2+G1) 

(IIB-10) 

L 

We  can  solve  equation  (IIB-10)  for  the  maximum  value  of 

S ,  the 

throughput.  We  obtain  Smax  =  .128  when  Gl  =  .71. 

L 

k 

To  include  the  effects  of  passive  acknowledgments 

we  use 

t 

t 

equations  (IIB-2  to  6)  to  obtain  for  arbitrary  length  chains, 

r 

S12  =  S 

f 

f* 

Sl'1+1  ^2,3^3,4--^i,i+l  '  1  “ . n  “ 

(IIB-11) 

► 

► 

Sn-l,n  =  Sn-2,n-l 

► 

► 

i 

L . 

•: 

J 
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Also, 

S12*  =  S23*  =  s 

s  *  =  S  (IIB-12) 

i,i+l  cf23...qi_1>i 

Returning  to  the  example  for  n=4  we  use  equations  (IIB-11)  and 
(IIB-3)  in  (IIB-1)  to  obtain 

S  =  G12P1 ,2,3)  =  q23G23P(l/2,3,4)  =  q23G21P(l  ,2,3) 

=  G12P(1,2,3)  =  G23  =  G21P1,2,3,4)  (IIB-13) 


or  G2  =  G1(l+G1)+G1(l+G1)2  =  G1(l+G1)(2+G1) 

Thus 

Gx  Gx 

S  =  1+Gx+2G2  =  (1+G1)+2G1(1+G1)(2+G1)  (IIB-14) 

The  throughput  here  is  Smax  =  .098  when  Gx  =  .37,  a  reduction  of 
23%. 

In  Table  IIB-1  we  compare  the  effects  of  passive  acknowledg¬ 
ments  for  various  length  chains.  We  see  that  as  the  chain  becomes 
long  the  throughput  vanishes.  Also  shown  in  that  table  is  the 

Table  IIB-1.  Effect  of  Passive  Acknowledgments 


Length 
of  Chain 

Perfect 

Acknowledgments 

Throughput 

Duplicates 

Detected 

Duplicates 
Not  Detected 

4 

.128 

.106 

.098 

5 

.111 

.083 

.069 

6 

.102 

.072 

.053 

7 

.097 

.066 

.044 

8 

.094 

.063 

.033 

9 

.092 

.061 

.034 

10 

.091 

.060 

.031 

00 

.086 

.057 

0 
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throughput  for  chains  when  duplicate  packets  are  detected  and  not 
transmitted  to  the  next  node.  We  analyze  this  case  now. 

If  duplicate  packets  are  detected  and  not  transmitted  further, 
then  S.  j*  is  just  the  desired  traffic  to  be  sent  from  i  to  j.  As 
before  we  find  S.  .*  from  the  offered  traffic  and  the  routing.  For 
example,  in  a  chain  S?  .+1  =  Sj*+^  .  =  S  for  i  =l,...,n-l.  Tnus 
equations  (JIB-1  to  3)  combine  to  give 


[P(N- ,N-]; 
c  *  =  G  _ - _ i— J— 

bi,j  ui,j  mr 


(IIB-16) 


except  at  terminal  nodes  where  no  further  transmission  is  necessary 

and  q..  =  1.  This  can  be  solved  as  before  to  yield  the  G-  •  as  a 
*J  *  /  J 

function  of  the  S-  .*  and  then  the  maximum  S.  .*  or  throughput  can 
i.J  M 

be  found.  Details  of  the  procedure  have  already  been  reported. 

We  present  the  analysis  for  the  four  node  chain  analyzed  above. 
Here  S.  ^  j  =  S  for  i  =  1,2,3.  Again  we  use  symmetry  to 

obtain 


(1+Gj) 

s  ~  Gj.  - n - G2 3 


1  _  ^  1+Gl 
ri+GjD  "  C'21  D~ 


(IIB-17) 


D  =  (1+Gi)2  +  ZGafl+GO 


as  before. 


But  G2  =  Gi  +  Gid+GO2,  so 


(1+Gi  ;+2Gj  [l+fl+Gi 


Solving  we  get  =  .106  when  G,  =  .42 


(IIB-18) 


IIB.7 


The  throughput  for  various  length  chains  are  given  in 
Table  IIB-1.  We  see  here  that  there  is  a  34%  reduction  in  throughput 
for  long  chains  when  the  effects  of  passive  acknowledgments  are 
included  and  duplicate  packets  are  detected  and  not  transmitted 
further. 

We  can  obtain  an  asymptotic  solution  for  the  infinitely  long 
chain.  First  assume  perfect  acknowledgments.  We  have  already 
discussed  this  in  a  previous  report,  but  present  a  different  analysis 
here.  We  assume  S.  .+^  =  S.+1  .  =  S  and  G.  -+^  =  G.+1  •  =  ^  G.  Then 


S  =  ^  G  P[i,i+l,i+2,i+3] 

_  1  r  SPfl _ ,i-l  ,i+4, . . .  ,n] 

2  u  SP(.l , . . .  ,nj 


(IIB-19) 


Here  SP[  ]  is  the  sum  of  products  notation  previously  introduced. 
We  have 

SP[l,...,i-l,i+4 . n]  =  SP[1 . i-l]SP[i+4 . n] 

and 

SP[1 . n]=GSP[i . i-l]SP[i+3 . n]+SP[l . i]SP[i+2 . n] 

Thus  S  =  G/2X  where 


SP 

fi+3,.. 

•  ,n' 

+  SPfl . i] 

SP[i+2, . . 

-nl 

SP 

(i+4,.. 

•  ,n 

\  [1 . i-1  J 

SP[i+4, . . 

•  ,n] 

(IIB-20) 


Furthermore, 

SP[i+2 . n]  =  GSP[i+4 . n]  +  SP[i+3 . n] . 


Therefore  if  we  let 


Q  = 


SP[1 . i] 

SPfl . i-1] 


SP[i . n] 

SP[i+l - n] 


(IIB-21) 


[- 

.  < 

1 13 . 8 

■  ■  '  '  ■  n 

then  asymptotically  for  large  i  and  n,  we  obtain 

( 

X 

=  GQ  +  Q(G+Q)  =  2GQ  +  Q2 

• 

or 

S 

G 

2(Q“+2GQ) 

(IIB-22) 

, 

We  can  solve  for  Q  by  noting  that 

SP[1 . i]  =  GSP[1 _ i-2]  +  SP[1 . i-1]  • 

i 

Thus 

• 

>. 

Q 

< 

< 

or 

G 

=  Q(Q-l) 

(IIB-23) 

; 

Note  that  as  G  goes  from  0  to  #,  Q  goes  from  1  to  ». 

We  finally 

obtain 

\9i 

k 

S 

_  Q-l 

2QC2Q-1) 

(IIB-24) 

l 

This  has  a 

maximum  of  Smax  =  .086  at  Q  =  1+1//2  or  G  = 

(l+VD/2. 

* 

If  the 

effect  of  passive  acknowledgments  is  included  in  the 

•  ■  : 

preceding  analysis  we  get 

S 

_lr  Pfi,i+l,i+2,i+312 

2  u  P[i+1  ,i+2,i+3] 

(IIB-25) 

But 

• 

■M 

* 

P 

[i,i+l,i+2,i+31  _  SPfl . i-1]  _  1 

(IIB-26) 

• 

P 

[i+1  ,i+2 ,i+3 j  SP[1, . . .  ,iJ  Q  ' 

; 

Thus 

i 

S 

_  Q-l 

2Q7(2Q-1) 

(IIB-27) 

* 

•: 

This  has  a  maximum  SmQX  =  .057  when  Q  =  1.4  and  G  =  .56 

* 

.  i 

■  9: 

_ ^  -  - 

. M 
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3.  Star  Networks 

We  next  analyze  the  effect  of  passive  acknowledgments  on 
throughput  in  star  networks.  We  will  consider  two  topologies  -  one 
where  all  legs  are  unconnected  from  each  other,  the  second  where  the 
nodes  closest  to  the  center  are  fully  connected.  We  have  studied 
these  networks  in  our  previous  work  and  assumed  there  that  all 
passive  acknowledgments  were  heard.  We  assume  further  here  that 
duplicate  packets  are  detected  and  not  retransmitted.  See 
Figure  IIB-2  for  the  configuration  of  the  star  netwrrk  with  zero 
connectivity,  L  legs,  and  K  nodes  in  each  leg.  For  full  connectivity 
the  nodes  0,1,K+1,2K+1, . . .  ,(L-1)K+1  are  fully  connected. 

The  equations  used  in  the  previous  section  are  still  valid  here. 
For  j,  a  neighbor  of  i, 

si.i  =Gi,j  p<tVV  (IIB'1) 

V  ■  «i.j  si,i  <IIB-2> 

qi,j  =  PCNj.NjJ/PCNj)  (IIB-3) 

When  duplicate  packets  are  detected,  all  S.  .*  =  S,  the  throughput  in 

*  /  J 

each  leg.  The  end  conditions  are 

q2K+L-l,(£+l)K  =  1  '  £  =  °' - L"1  (IIB-28) 

since  transmissions  to  the  end  nodes  are  always  heard  and  passive 
acknowledgments  are  not  sent.  This  can  be  found  from  equation 
(IIB-3).  We  further  assume  that 

q£K+i,0  =  l  ,  I  =  0 . L-l  (IIB-29) 

Transmissions  to  the  center  node  are  not  always  heard.  Retransmis¬ 
sions  due  to  collisions  are  included  in  our  analysis.  However  since 
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the  center  node  does  not  retransmit,  passive  acknowledgments  are  not 
available.  We  assume  a  perfect  end-to-end  acknowledgment  is 
operating. 

Equations  (IIB-1  to  3)  have  been  solved  with  perfect  passive 
acknowledgments  and  with  imperfect  passive  acknowledgments  and 
duplicate  packet  detection.  Two  topologies  have  been  studied  -  zero 
connectivity  and  full  connectivity  around  the  center.  The  resultant 
maximum  throughputs  per  leg  are  shown  in  Table  IIB-2  and  3.  The 
maximum  throughput  from  the  center  node  is  given  in  Table  IIB-4. 
We  see  from  these  results  that  there  is  again  a  34%  reduction  in 
throughput  for  long  legs  and  zero  connectivity.  For  full  connectivity 
we  find  a  43  to  46%  reduction  in  throughput  due  to  the  failure  to  hear 
passive  acknowledgments . 

From  our  previous  work  we  found  that  except  for  short  legs,  the 
fully  connected  topology  is  better  when  passive  acknowledgments  are 
not  considered.  This  is  because  inward  transmissions  from  neighbors 
of  the  center  are  not  collided  with.  With  zero  connectivity  passive 
acknowledgments  for  transmissions  from  the  first  ring  to  the  second 
are  interfered  with  only  by  the  center  node.  With  full  connectivity 
all  other  first  ring  nodes  add  to  this  interference.  Thus  the  first 
ring  will  have  to  retransmit  more.  In  Table  IIB-4  we  compare  the 
maximum  throughputs  from  the  center  node  for  full  and  zero  con¬ 
nectivity  and  see  that  indeed  zero  connectivity  is  better. 


Table  IIB-2:  Effects  of  Passive  Acknowledgements  on  Star 

Networks.  Zero  Connectivity  (Number  of  Legs, 
L=5)  -  Duplicates  Detected 


Maximum  throughput  per  Leg,  S 

Number  of  Hops  - 

er  Leg,  K  Perfect  Acknowledgments  Imperfect  Acknowledgments 
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1 

A 

> 


Table  IIB-4 

:  Effect  of 

Maximum 

Passive  Acknowledgments  on 
Throughput  from  the  Center. 

Star  Networks. 

4 

►  .  •; 

Number  of 
Legs,  L 

Number  of 
Hops  per 
leg,  K 

Connectivity 

Maximum  throughput  per  Leg,  Smax 

Perfect  Imperfect 

Acknowledgments  Acknowledgments 

i 

4 

>  • 

5 

2 

Zero 

.225 

.165 

>3 

.220 

.145 

—1 

6 

3 

FuH 

.240 

.140 

►  •: 

>4 

.240 

.138 

- 

9 

3 

Fun 

.257 

.140 

j 

>4 

.252 

.139 

>  mi 

4.  A  new  protocol  for  acknowledgments 

In  previous  sections  we  studied  the  effects  of  passive  acknowl¬ 
edgments  on  the  maximum  obtainable  throughput  for  various  topologies. 
We  found  that  the  throughput  decreased  by  34  to  46%.  Disadvantages 
of  the  proto -.ol  studied  are: 

1)  each  node  has  only  one  chance  to  hear  acknowledgments 
from  its  neighbors, 

2)  failure  to  hear  passive  acknowledgments  causes  increased 
channel  traffic  since  packets  must  be  repeated, 

3)  the  next  node  must  recognize  duplicate  packets  or  perfor¬ 
mance  will  be  drastically  worse,  and 

4)  a  node  must  wait  a  time-out  period  before  retransmitting, 
thus  incurring  further  delays. 


>  ■  mi 


>  mi 


l  9i 


1 


i 

i 
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We  propose  the  following  protocol.  A  node  in  general  has  many 
packets  to  transmit,  in  addition  to  repeating  a  packet  just  received. 
These  packets  are  intended  for  all  its  neighbors.  The  old  protocol 
requires  one  to  wait  until  that  particular  packet  is  repeated.  Here 
we  enlarge  the  header  of  every  packet  so  that  it  includes  acknowl¬ 
edgment  information  for  all  recently  received  packets.  Thus  a 
"passive"  acknowledgment  can  be  received  on  the  next  packet  trans¬ 
mitted.  The  acknowledgment  for  a  particular  packet  can  be  included 
in  the  next  several  packets,  so  there  are  several  attempts  to  hear  the 
acknowledgment.  At  the  expense  of  increased  header  and  more  logic 
in  the  nodes,  the  effect  of  passive  acknowledgments  can  be  made 
negligible. 

We  describe  the  protocol  more  carefully.  A  node  i  sends  a 
packet  to  one  of  its  neighbors,  node  j.  This  packet,  if  heard,  is  to 
be  relayed  to  a  neighbor  of  j,  say  k.  Node  j  has  other  neighbors. 
The  next  m  transmissions,  after  receipt  of  the  packet  in  question, 
from  node  j  to  any  of  its  neighbors  includes  in  the  header  a  specific 
acknowledgment  for  the  packet.  Node  i  is  listening  and  has  m 
chances  to  hear  the  acknowledgment.  Node  i  sets  a  time-out  period, 
which  is  the  expected  time  for  node  j  to  transmit  m  times.  Note  that 
the  acknowledgment  is  in  general  received  (or  the  time  out  period 
reached)  before  the  packet  is  actually  relayed.  A  prelimirtary 
analysis  of  this  protocol  indicates  that  for  m  =  3  or  4  the  effect  of 
passive  acknowledgments  can  be  reduced  drastically  for  many  topol¬ 
ogies.  Full  details  will  be  presented  in  the  next  report. 
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IIC.  Switching  For  a  Large  Number  of  Users  With  Short 

Packets  and  a  Tight  Time  Delay  Constraint. 

1.  Introduction 

One  of  the  main  areas  of  our  research  is  to  investigate  new 
communications  techniques  within  networks.  The  motivation  for  this  it 
two-fold.  First,  to  take  advantage  of  the  vast  increases  in  processing 
power  and  storage  now  available  in  networks.  Second,  packet  switch¬ 
ing  has  been  found  to  perform  well  for  certain  classes  of  problems 
but  not  for  all.  Other  techniques  are  better  suitable  for  file  transfers, 
facsimile,  integrated  data  and  voice,  very  short  packets,  telemetered 
data,  etc.  For  very  long  messages  some  form  of  circuit  switching  is 
preferable.  For  non-bursty  traffic  some  form  of  dedicated  assignment 
is  preferable.  Finally,  for  very  short  messages  the  overhead  in 
packet  switching  becomes  burdensome. 

We  concentrate  here  on  the  latter  problem.  In  this  report  we 
will  analyze  the  performance  of  packet  switching  techniques  and 
develop  optimum  packetizing  strategies  for  this  problem.  In  subse¬ 
quent  reports  we  will  develop  and  analyze  more  efficient  communica¬ 
tion  schemes. 

The  problem  we  consider  is  that  of  a  large  number  of  users 
communicating  from  remote  terminals  to  a  host  computer.  They  input 
characters ,  one  at  a  time ,  in  a  bursty  fashion  -  short  ’  bursts  of 
characters  with  varying  times  between  characters,  followed  by  long 
pauses  between  bursts.  The  echoing  of  characters  is  performed  by 
the  host  computer.  The  round  trip  delay  from  the  time  a  character  is 
typed  until  it  is  received  must  be  short  or  further  input  is  impeded. 


One  solution  (a  common  one)  is  to  form  a  packet  for  each  charac¬ 
ter.  Each  packet  must  contain  a  header  which  includes  synchroniza¬ 
tion  information,  addresses  of  the  large  number  of  users,  error 
checking,  etc.  Thus  the  overhead  is  excessive.  Longer  packets 
reduce  the  overhead  but  at  the  expense  of  packetization  delay.  We 
will  develop  and  analyze  optimum  packet  switching  strategies  as  a 
prelude  to  investigating  other  switching  techniques. 

In  our  model  we  consider  ail  users  connected  to  a  concentrator 
which  is  then  connected  to  the  host  computer.  Except  for  the  echo¬ 
ing  we  ignore  the  return  messages  from  the  host  to  the  users.  The 
total  delay  experienced  includes  packetizing  delay,  queuing  delay  in 
the  concentrator,  and  transmission  time  over  the  link  to  the  host. 
Other  delays,  such  as  host  turnaround  time,  are  considered  either 
negligible  or  invariant  to  the  schemes  we  compare. 

Our  goal  is  to  find  the  maximum  number  of  users  that  can  com¬ 
municate  with  an  allowable  maximum  delay.  To  that  end  we  analyze 
and  compare  various  packetization  strategies  below.  We  also  dev- lop 
a  packetization  strategy  which  is  optimum. 

2.  Model  Descriptions 

First  consider  one  user.  It  generates  characters  in  a  bursty 
fashion  with  long  idle  periods.  We  consider  time  slotted  in  A  second 
periods.  (A  is  the  minimum  time  between  characters,  say  1/4  second.) 
When  the  user  is  generating  characters  we  assume  it  produces  a 
character  in  a  slot  with  probability  p.  With  probability  1-p  there  is 
no  character  in  that  slot.  Successive  slots  are  independent.  Thus 
the  rate  of  generating  characters,  when  busy,  is  p/A  characters  per 
second.  (For  p  =  1/4  and  A  =  1/4  this  becomes  one  character  per 
second.)  The  average  time  between  characters  is  A/p  seconds. 
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The  characters  are  collected  .n  a  buffer  until  a  packet  is 
formed.  Then  a  header  is  added  to  the  packet,  which  is  sent  over 
an  error-free  and  synchronous  channel  with  capacity  C  characters  per 
slot.  The  packet  is  received  by  the  host,  which  then  echoes  the 
characters,  i.e.,  sends  tnem  back  to  the  user.  They  then  appear  on 
the  user's  terminal  as  confirmation  of  proper  receipt. 

Vve  ignore  the  return  message  from  host  to  user  in  this  analysis, 
and  consider  only  the  inbound  packets  and  echoed  characters.  We 
consider  only  the  delays  in  packetization  and  queuing  (at  the  buffer) 
and  in  transmission  over  the  line  to  the  host.  Other  delays,  such  as 
host  turnaround  time,  are  either  negligible  or  invariant  to  the  schemes 
we  compare.  We  also  ignore  the  effect  of  users  switching  between 
active  and  idle  states. 

Our  model  is  shown  in  Figure  IIC-1. 


character  packets  characters/ 

stream  formed  slot 


Figure  IIC-1.  A  single  user  model. 


We  are  interested  in  finding  the  number  of  users,  such  as  the 
one  described  above,  that  can  be  handled  by  the  system  with  tolerable 
delay.  Assume  that  an  any  time  (or  for  a  short  period  of  many  slots) 
M  of  a  larger  population  of  identical  users  are  actively  generating 
characters.  We  want  to  find  the  maximum  tolerable  value  for  M.  The 
size  of  the  larger  population  influences  only  the  size  of  that  part  of 
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the  header  needed  for  addressing.  This  is  not  sensitive  and  will  not 
be  considered  further  here.  _ 

We  model  each  of  the  users  the  same  way  we  modeled  the  single 
user  above.  The  M  users  are  generating  characters  independently  of 
each  other.  The  characters  are  collected  in  the  buffer  and  formed  .  - 

into  separate  packets  for  each  user.  They  are  then  queued  and 
transmitted. 

This  model  is  shown  in  Figure  IIC-2.  ^ 


Figure  IIC-2.  A  Multi-User  Model. 


3.  Analysis  of  Fixed  Length  Packet  Switching 

We  first  analyze  a  fixed  length  packet  switching  scheme  for  a 
single  user.  Here  the  buffer  waits  until  y  characters  are  received 
and  then  forms  a  packet  by  added  h  characters  as  header.  Thus  the 
packet  length  is  x  =  y+h  characters.  The  delay  we  are  concerned 
with  is  given  by 


D  =  PD  +  W  +  S  (IIC-1) 

where  D  is  the  total  average  delay  per  character, 

PD  is  the  average  packetization  delay  per  character, 

W  is  the  average  waiting  time  for  the  packet,  and 
S  is  the  service  time  or  transmission  time  for  the  packet. 
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All  times  are  in  slots. 

For  fixed  length  packets  of  length  x,  S  =  x/C  =  (y+h)/C.  To 
find  the  average  packetization  delay  note  that  the  average  time 
between  generation  of  characters  is  1/p  slots.  The  packet  begins  to 
be  formed  upon  generation  of  the  first  character  and  is  considered 
formed  when  the  y^  character  is  in  the  buffer.  Thus  the  first 
character  must  wait  on  the  average  (y-l)/p  slots,  the  last  character 
does  not  wait  at  all,  and  the  k111  character  waits  (y-k)/p  slots.  The 
average  packetization  delay  per  character  is 


or 


PD  = 


PD  = 


I  /  ZdS 
*k=l  p 


Yil 

2p 


(IIC-2) 


Other  parameters  of  interest  are 

a  =  average  time  between  packets  =  y/p 
A.  =  packet  arrival  rate  =  1/a  =  p/y 
p  =  line  utilization  = 


To  find  the  average  waiting  time  we  must  first  find  the  distribu¬ 
tion  of  the  time  between  packets.  The  time  between  characters  is 
geometric  and  has  moment  generating  function  pz/(l-qz),  where 
q=l-p.  The  y  successive  times  between  packets  are  independent,  so 

A(z)  =  moment  generating  function  of  time  between 
packets 


1-qz 


y 


(IIC-3) 
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The  distribution  of  the  time  between  packets  is  given  by 

P[a=£]  =  (y_\)  q£"y  pY  ,  £  >  y.  (IIC-4) 

00  2 

This  can  be  found  by  using  A(z)  =  I  P[a=£]z  and  expanding  equation 

£=y 

(IIC-3)  or  by  noting  that  the  £  slots  after  a  packet  is  formed  must  be 
filled  by  y-1  characters  in  the  first  £-1  slots.  The  last  (£t^1)  slot 
must  contain  the  last  character. 

Consider  the  problem  of  packets  arriving  to  a  queue  with  inter¬ 
arrival  times  given  by  equations  (IIC-3)  or  (IIC-4).  The  service 
length  is  a  constant  given  by  S  =  (y+h)/C  slots.  This  is  a  G/D/l 
queue.  The  generating  function  for  the  waiting  time  distribution  can 
be  found  as 


W(z)  = 


N(z) 


Kz-q)y-pyzX] 


(IIC-5) 


x  y-i 

where  N(z)  =  q"  2 

i=0 


W(l)(q) 

i! 


y-!-!  x  1 
2  (X  )^  R  (z) 

j=0  J  q  J 


R.(z)  =  Pi+^(z-q)y  -  Py(z-q)i+^ 


x  =  y+h  and  q  =  1-p 


The  average  waiting  time  is  given  by 


W 


,y-2. 


=  w(1)(i)  -N"(1)~P:r  ‘Jy(y-D+Pyx(x-i 

2py-1(y-px) 


-l) 


(IIC-6) 


where  N"(l)  = 


d2N(z) 

dz2 


Y-1  wO), 


y-i-i 


Z=1 


i=0 


1  ’ 


i=i 


1  rJ 


[y(y-i)  -  (i+i)  (i+j-D]  1 


We  can  solve  for  W  by  the  following  method.  First  find  :h« 
(y-1)  roots  of  (z-q)y-pyzx  =  0  (from  equation  IIC-5;  which  are  insult 
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the  unit  circle,  i.e.,  |z|  <  1.  For  these  values  of  z,  the  numerator 
N(z)  must  vanish.  Thus  N(z)  =  0  for  these  y-1  values  of  z.  This 
gives  y-1  linear  simultaneous  equations  for  the  y  values  of  W^(q). 
The  last  equation  is  obtained  by  noting  that  z  =  1  is  also  a  root, 
Rj(l )  =  0,  so  N(z)  =  0  also.  But  W(l)  =  1.  L’Hospital's  rule  can  be 
used  to  find  this  last  relation. 

The  average  (inbound)  character  delay  can  now  be  found  from 
equation  (IIC-1). 

4.  Variable  Length  Packet  Switching  Techniques 

In  this  section  we  present  three  different  variable  length  packet 
switching  techniques.  In  all  cases  we  use  the  same  user  model  as 
presented  above.  In  the  next  sections  we  will  analyze  the  perfor¬ 
mance  of  these  techniques  and  compare  them  with  fixed  length  packet 
switching . 

Variable  Length  Packet  Switching  Technique  I  (VLPS  I). 

In  this  technique  a  packet  is  generated  every  r  slots.  However 
if  any  r  slot  interval  contains  no  characters,  then  the  next  r  slot 
interval  starts  when  the  next  character  is  produced. 


» 


I 


I 


»  • 


*  -mi 


•  -•< 


Variable  Length  Packet  Switching  Technique  II  (VLSP  II).  I 

Here  a  packet  is  generated  every  r  slot.  If  any  r  slot  interval 
is  empty  then  it  is  skipped  and  the  next  r  slot  interval  is  considered. 

The  two  schemes  differ  in  the  following  manner.  In  VLSP  II  p 

time  is  blocked  into  r  slot  units.  A  packet  is  sent  at  the  end  of  the 
block  if  at  least  one  character  has  been  produced.  In  VLSP  I,  when 
a  block  has  no  characters,  the  next  block  is  not  formed  until  a  new  p 

character  is  produced.  A  third  scheme,  in  which  all  blocks  start 
when  the  first  character  is  produced,  is  not  considered  here. 


r 


Modified  Fixed  Length  Packet  Switching  (MFLPS) 

A  packet  is  produced  upon  the  arrival  of  either  n  or  n+1  char¬ 
acters.  The  choice  is  made  randomly.  Thus,  let  1-Q  be  the  prob¬ 
ability  that  the  packet  contains  n  characters.  With  probability  Q,  the 

s  t 

packet  will  be  formed  when  the  next  (n+1 — )  character  arrives.  The 
packet,  of  either  length,  is  created  immediately  upon  the  arrival  of 
the  last  character.  We  have  proven  that  this  scheme  is  optimum  in 
the  sense  that  for  a  fixed  average  packetization  delay  it  yields  the 
maximum  average  packet  length.  We  will  describe  this  further  in  the 
next  report.  We  will  describe  in  the  next  section  a  procedure  for 
finding  the  optimum  values  of  n  and  Q. 

The  fixed  length  packet  switching  technique  (FLPS)  was  de¬ 
scribed  and  analyzed  in  the  previous  section. 


5.  Analysis  of  Variable  Length  Packet  Switching  Techniques 
Variable  Length  Packet  Switching  Technique  I 

The  following  are  the  results  of  our  analysis.  Most  of  the 
parameters  have  already  been  defined  in  section  3. 

—  r  "I 

y  =  average  packet  length  =  (1+q  )rp  +  q 
a  =  average  time  between  packets  =  y/p. 

A  =  arrival  rate  of  packets  =  1/a. 


q  =  1-p 

_  V'.U.  ”1 

y2  =  second  moment  of  packet  length  =  ap(pr+q)  +  q  p(r-l) 
x  =  y  +  h 
x  -  y  +  h 


72  = 


=  y2  +  2hy  +  h2 


PD  =  average  packetization  delay 

=  A  {  0>l)(r-2)(l+q_).  +  (r-l)(l  +  )  } 

P 


+ 


( 1 1 C  -  7  ) 
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S  is  the  average  service  time  per  character.  Although  each 
character  in  a  packet  must  wait  for  the  entire  packet  to  be  received 
we  consider  that  service  time  to  be  apportioned  over  ail  characters  in 
the  packet.  A  packet  with  y  characters  has  a  service  time  of  (y+h)/c. 
Thus 

S  =  HnEEI  =  if  ♦  fc  .  (IIC-?a) 

y  y 

We  consider  M  identical  users.  We  approximate  their  waiting  time 
at  the  buffer  by  assuming  it  acts  as  an  M/G/l  queue.  Thus  the  total 
average  delay  per  character  is  given  by 


15  =  PD  +  s  ^TFp)  (IIC"8)  * 

where  p  =  utilization  =  <  1  . 

Equation  (IIC-8)  can  be  written  so  as  to  give  M  as  a  function  of 
D.  We  will  find  the  value  of  the  parameter  r  to  maximize  M  for  a 
fixed  D. 

Variable  Length  Packet  Switching  Technique  II. 

Using  the  same  definitions  as  above,  the  results  of  the  analysis 
of  this  scheme  are  given  by: 

a  =  —  =  1/A 
l-qr 


Y2  ~  (rp+q) 

1-q 


PD  = 


r-1 

2 


(IIC-9) 
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We  can  again  use  Equation  IIC-8  to  find  D  or  to  find  the  maximum 
M  for  a  fixed  value  of  D. 

Modified  Fixed  Length  Packet  Switching 

Recall  that  the  packet  length  is  either  n,  with  probability  1-Q, 
or  n+1  with  probability  Q.  The  packets  are  created  upon  receipt  of 
the  last  character.  The  values  of  n  and  Q  are  to  be  determined. 
(We  have  already  mentioned  the  optimality  of  this  technique.)  The 
results  of  the  analysis  are: 


y  =  n+Q 


(IIC. 10) 

T2  =  (n+Q)2  +  Q(l-Q) 

_  Mp(n+Q+h) 
p  C(n+Q) 

From  these  results  we  can  again  get  an  expression  for  D,  and 
for  M.  We  can  then  find  the  values  of  the  parameters  y,r,n  and  Q  in 
the  various  schemes  to  maximize  M  for  a  given  value  of  D.  This  is 
done  in  the  next  section. 

6.  Comparison  of  Different  Techniques 

Each  of  the  four  schemes  described  above  have  been  analyzed. 
For  each  scheme  we  found  the  optimum  value  of  their  parameters 
(y,r,n,Q  respectively)  so  as  to  maximize  M  for  a  given  value  of  D. 
Thus  we  seek  the  maximum  number  of  users  that  can  be  handled  while 
providing  an  allowable  total  character  delay. 

We  have  considered  an  example  where  p  =  1/4,  C  =  300  char¬ 
acters  per  slot,  and  a  slot  is  1/4  second  long.  Thus  the  line  speed 
is  1200  characters  per  second.  One  source  produces  on  the  average 
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one  character  per  second  with  a  peak  of  4  per  second.  Without 
overhead  and  ignoring  delay  at  most  1200  users  can  be  handled. 
However  both  overhead  and  delay  constraints  reduce  this.  We  have 
considered  header  sizes  of  1  or  2  characters  per  packet  and  a  maxi¬ 
mum  delay  of  D  =  4  slots  or  1  second.  The  results  are  shown  in 
Table  IIC-1. 


Table  IIC-1  Comparison  of  Different  Packet  Switching  Techniques 


Optimum  Value 

of  Design 

Technique 

Parameter 

a 

Y 

PD 

h 

M 

Fixed 

Length  PS 

*< 

ii 

8 

2 

2 

1 

797 

800 

2 

597 

600 

Variable 

Length  PSI 

r=8 

9.1 

2.27 

3.61 

1 

821 

833 

2 

626 

638 

Variable 

Length  PSII 

r=8 

8.89 

2.22 

3.5 

1 

817 

827 

2 

622 

631 

Modified  Fixed 

Length  PS 

n=2,Q=0.77 

11.08 

2.77 

3.67 

1 

864 

881 

n=2,Q=0.79 

11.16 

2.79 

3.70 

2 

679 

698 

We  can  see  from  the  results  that  the  modified  fixed  length 
packets  switching  scheme  performs  best.  Indeed  we  have  proven  that 
this  scheme  is  optimum  when  considering  only  packetization  delay  and 
service  time.  The  variable  length  schemes  perform  about  the  same 
and  better  than  the  fixed  length  scheme.  The  value  M*  in  the  tables 
is  that  for  which  p  =  1  and  serves  as  an  upper  bound  on  M.  As  M 
approaches  M*  the  waiting  time  dominates  and  becomes  large. 


lie. 12 


7.  Conclusions 

We  have  analyzed  four  different  packet  switching  techniques  in 
an  environment  where  packet  switching  is  an  inefficient  technique. 
We  considered  many  users  generating  individual  characters  and 
requiring  a  fast  response.  Small  packets  waste  overhead,  while  large 
packets  consume  the  response  time.  The  optimization  is  very  tight. 
We  have  found  a  modified  fixed  length  packet  switching  technique 
which  has  some  optimal  properties  and  performs  the  best  of  those 
schemes  studied. 

Subsequent  work  will  further  explore  the  optimality  of  the 
modified  fixed  length  packet  switching  technique,  and  explore  the 
performance  of  packet  switching  techniques  over  a  wide  variety  of 
situations  and  design  parameters  --  source  models,  time  delay  con¬ 
straints,  header  sizes,  etc. 

We  intend  to  focus  attention  on  other,  more  sophisticated,  techniques 
to  transmit  the  characters  from  many  users.  Thus  we  seek  alter¬ 
natives  to  packet  switching  in  these  environments. 


IID.  Probabalistic  Analysis  of  Algorithm  Performance 

The  problem  of  finding  an  optimal  set  of  repeater  locations 
covering  a  given  set  of  terminal  sites,  or  set  of  potential  repeater 
locations  and  a  covering  matrix  specifying  which  terminal  sites  can  be 
covered  by  each  repeater  location,  is  an  important  problem  in  the 
design  of  multihop  packet  radio  networks.  The  problem  in  an  instance 
of  the  classic  set  covering  problem  which  has  many  other  applications 
as  well.  The  set  covering  problem  is  known  to  be  NP-complete  and 
as  such  it  is  unlikely  that  an  algorithm  will  be  found  which  can 
guarantee  an  optimal  solution  and  also  guarantee  a  reasonable  running 
time,  that  is  a  runtime  which  grows  polynomially  in  the  number  of 
terminals  and  repeaters. 

The  theory  of  nondeterministic  polynomial  completeness  (NP- 
completeness)  basically  states  that  there  is  a  large  class  of  problems, 
which  includes  almost  all  optimization  problems  of  interest  in  the  area 
of  network  design,  which  have  the  property  that  a  solution  to  one 
could  be  used  to  obtain  a  solution  to  all  the  others  in  a  reasonable 
amount  of  time.  It  has  been  shown  that  one  problem,  the  Satisfiability 
Problem,  has  the  property  that  all  problems  in  this  class  can  be 
transformed  into  an  instance  of  it.  Thus,  an  algorithm  which  could 
obtain  an  optimal  solution  to  the  Satisfiability  Problem,  could  be  used 
to  solve  all  the  others.  Techniques  have  been  developed  for  proving 
that  many  other  problems  are  NP-complete.  Proofs  have  been  given 
for  the  NP-completeness  of  over  200  other  well-known  problems. 

In  our  previous  work  we  reported  on  results  obtained  using  a 
new  set  covering  algorithm.  The  algorithm  obtained  optimal  solutions 
for  moderate  sized  problems  with  up  to  300  node  networks  and  cover- 
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age  related  to  distance.  That  is,  if  the  average  matrix  reflected  the 
fact  that  a  repeater  was  far  more  likely  to  cover  a  nearby  terminal 
than  one  farther  away,  then  the  algorithm  worked  very  well,  con¬ 
verging  to  an  optimal  solution  very  quickly. 

If  on  the  other  hand,  the  coverage  matrix  was  random,  the 
algorithm  ran  much  more  slowly  and  only  problems  with  up  to  50 
terminals  could  be  treated  efficiently.  In  examining  why  the 
algorithm  had  difficulty  treating  problems  with  random  data,  we 
observed  that  the  difficulty  stemmed  from  the  large  number  of  optimal 
and  near  optimal  solutions  available.  The  algorithm  had  no  difficulty 
finding  an  optimal  solution  but,  rather,  had  difficulty  in  verifying  the 
solution  was  optimal  in  a  reasonable  amount  of  time  due  to  the  presence 
of  a  large  number  of  alternate  solutions  of  comparable  quality.  This 
led  us  to  conjecture  that  it  might  be  possible  to  develop  algorithms 
which  had  a  reasonable  running  time  and  which  could  find  optimal 
time  or  near-optimal  solutions  with  high  probability.  It  is  known  that 
asymptotically,  as  the  size  of  the  problem  becomes  infinite,  there 
exists  algorithms  which  give  optimal  results  almost  always.  These 
results  are  an  outgrowth  of  the  same  symmetry  which  makes  it  difficult 
to  verify  the  optimality  of  the  solutions  produced  by  our  set  covering 
algorithm. 

These  results  are  encouraging  but  of  no  direct  use  in  solving 
network  design  problems  since  real  problems  are  of  finite  size  and 
nothing  was  said  about  how  fast  the  algorithms  converge  probablis- 
tically  to  an  optimum;  i.e.,  what  the  probability  is  of  their  finding 
the  optimum,  or  a  solution  within  some  bound  of  the  optimum.  We 
thus  sought  to  investigate  what  could  be  said  about  the  probabalistic 
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..  t^}  of  N  terminals,  a  set 
rM)  of  M  potential  repeater  sites,  and  an  N  x  M 


performance  of  known  heuristics  for  problems  of  moderate  size.  We 
chose  the  set  covering  problem  mentioned  above  as  the  first  problem 
to  investigate.  The  probabalistic  model  of  this  problem  follows. 

We  are  given  a  set  T={t1,  t2, 

R={rx,  r2,  .. 

covering  matrix,  C,  where  c.j  is  1  if  t.  can  be  covered  by  r.  We 
seek  a  subset  of  R  containing  as  few  repeaters  as  possible  and  cover¬ 
ing  all  terminals.  The  elements  of  C  are  chosen  randomly  and  inde¬ 
pendently.  Specifically,  we  assume  Pr{c.j=l}=p  and  Pr{C.j=0}=l-p=q. 
We  define  Rj,  the  set  of  terminals  covered  by  r^,  by: 

Rj  =  JtjICy  =  I) 

Similarly,  we  define  Tj,  the  set  of  repeaters  covering  t.,  by: 

Ti  =  fri'cij  *  n 

A  cover  of  the  terminals  is  then  defined  as  a  subset  SCR  satisfying: 


U  R.  =  T 
r.eS  ] 

A  minimum  cover,  S*,  is  then  a  cover  containing  as  few  elements  as 
possible,  i.e.: 


U 

r^eS* 


and  |S*|  <  | S |  for  all  S  such  that  U  R.  =  T. 

r.eS  ] 

Note  that  the  minimal  cover  is  not  in  general  unique.  Indeed 
our  approach  rests  on  the  existence  of  a  large  number  of  minimal  and 
near  minimal  covers. 

In  order  to  analyze  the  probabalistic  behavior  of  a  heuristic  we 
must  do  two  things.  First,  for  a  given  set  of  values  of  N,  M,  and 
p,  we  must  find  the  probability  that  |S*1  =  K,  that  is,  the  probability 
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that  there  exists  a  minimal  covering  containing  exactly  K  repeaters. 
Let 

P(N,M,K)  =  Pr{|S*|  =  K}  for  given  N,M,  and  p. 


Next  we  must  find  PA^G(N,M,K),  defined  as  the  probability  that  the 
algorithm  will  find  a  solution  with  exactly  K  repeaters.  A  figure  of 
merit  for  evaluating  the  probabilistic  behavior  of  the  algorithm  is 
then 


FALG 


M  M 

I  KP.Tr(N,M,K)  -  I  KP(N,M,K) 
K=0  K=0 


M 

I  KP(N,M,K) 
K=0 


(IID.l) 


This  is  simply  the  relative  error  made  by  the  algorithm,  i.e., 
the  difference  between  the  avarage  number  of  repeaters  in  a  covering 
found  by  the  algorithm  and  the  minimum  number  required,  normalized 
by  the  minimum  number  required.  We  believe  that  FALG  0  for 

any  M,  and  p,  and  any  reasonable  algorithm  where  a  reasonable 
algorithm  is  defined  as  one  which  stops  when  it  has  a  cover.  In 
particular,  we  believe  F^q^O  for  an  algorithm  which  picks  repeaters 
randomly  until  a  covering  is  obtained. 

We  begin  by  evaluating  Q(N,M,K),  the  probability  that  a  ran¬ 
domly  chosen  set  of  K  repeaters  covers  all  N  terminals.  Notice  that 
this  is  not  the  same  as  the  probability  that  an  algorithm  picking 
repeaters  at  random  will  stop  after  K  repeaters  since  in  some  cases 
fewer  than  K  will  also  suffice  and  the  algorithm  will  find  them.  Also, 
Q(N,M,K)  is  a  cumulative  distribution,  that  is  it  relates  to  the  proba¬ 
bility  that  K  or  fewer  repeaters  are  required  for  a  covering  rather 
than  exactly  K.  We  do  have,  however,  that 
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K 

Q(N,M,K)  <  I  P(N,M,K)  (IID.2) 

L=0 

that  is,  the  probability  of  a  randomly  chosen  set  of  K  repeaters 
covering  all  N  terminals  is  less  than  or  equal  to  the  probability  of  the 
existence  of  a  covering  containing  K  or  fewer  terminals.  Thus  we 
have, 

P(N ,M,K)  >  Q(N,M,K)  -  Q(N,M,K-1),  K>1, 
that  is,  the  difference  between  two  successive  Q's  is  a  lower  bound 
on  the  corresponding  P. 

Q(N,M,K)  is  easy  to  evaluate  because  repeaters  picked  randomly 

are  picked  independently  of  one  another  and  our  random  model 

assumed  their  original  characteristics  are  independent.  The 

probability  of  a  single  terminal  not  being  covered  by  a  single  repeater 

is  q(=l-p).  The  probability  of  a  single  terminal  not  being  covered  by 

any  of  the  K  repeaters  is  q  .  The  probability  of  at  least  one  repeater 

of  K  covering  a  given  terminal  is  then  1-q  .  Finally,  the  probability 

K  N 

that  all  N  terminals  are  covered  by  the  K  repeaters  is  (1-q  )  .  We 
thus  have  that 

Q(N,M,K)  =  (1  -  qK)N  (IID.3) 

Note  that  this  is  not  a  function  of  M  except  in  the  trivial  sense 
that  M  must  be  no  smaller  than  K.  Figure  IID.l  shows  values  of 
Q(N,M,K)  for  p=.5.  Figure  IID.2  shows  values  of  Q(N,M,K)  - 
Q(N,M,K-1)  for  the  same  examples.  Both  show  a  tendency  for  solu¬ 
tions  to  cluster  about  a  narrow  range  of  K.  This  is  encouraging  in 
that  it  implies  that  probabalistically ,  simple  algorithms  should  do  well. 
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FIGURE  IID.l  CUMULATIVE  PROBABILITY  OF  A  COVERING  FOR 
RANDOM  ALGORITHM  (p=0.5). 
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As  yet,  we  do  not  have  an  exact  analytic  expression  tor 
P(N,M,K).  We  can  get  a  tighter  lower  bound  on  it  by  evaluating  the 
distribution  of  solutions  from  increasingly  sophisticated  algorithms.  A 
simple  minded  algorithm  would  be  to  pick  repeaters  randomly  until  one 
obtained  a  covering.  The  probability  of  obtaining  a  solution  with  K 
or  fewer  repeaters  using  this  algorithm  is  P1(N,M,K)  given  by: 

N 

Pi(N,M,K)  =  Z  f  piqN"1P1(N-i,  M-1,K-1)  (IID.4) 

i=0  1 

since  the  K—  repeater  will  cover  i  terminals  with  probability  (^piq^-* 
and  the  remaining  K-l  of  M-l  repeater  must  cover  the  remaining  N-i 
terminals.  Note  that  this  is,  again,  not  a  function  of  M  except  that 
M  must  be  no  smaller  than  K.  The  probability  of  obtaining  a  covering 
with  exactly  K  repeaters  is  given  by  P1(N/M,K)-P1(N,M,K-1). 

A  somewhat  more  sophisticated  algorithm  is  to  consider  repeaters 
in  a  random  order  and  to  select  a  repeater  only  if  it  improves  the 
chances  of  obtaining  a  covering  more  by  picking  it  than  by  not 
picking  it.  Letting  P2(N,M,K)  be  the  probability  of  obtaining  a 
covering  containing  K  or  fewer  repeaters  using  this  algorithm,  we 
have 

N 

P2(N,M,K)  =  Z  N  p1qN"1{max  (P2(N-i,  M-l, K-l),  P2(N,M-1,K))} 
i=0  1 

(IID.5) 

P2(N,M,K)  is  plotted  for  N=5  and  N=10  in  Figure  IID.3.  Along  with 
it,  the  probability  of  obtaining  a  covering  using  K  randomly  chosen 
repeaters  is  plotted.  Notice  that  this  algorithm's  performance  is  a 
substantial  improvement  over  randomly  picking  repeaters.  Notice 
also,  that  P2(N,M,K)  is  a  function  of  M.  This  algorithm  is  very 
simple  and  has  a  reasonable  running  time  even  for  very  large  N  and 
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M.  Our  next  objective  is  to  compare  it  with  an  upper  bound  on 
P(N,M,K)  to  see  if  it  already  has  converged  to  near  optimal 
performance  probabalistically . 
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